Propagation inhibition and wave localization in a 2D random liquid medium 
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Acoustic propagation and scattering in water containing many parallel air-filled cylinders is studied. 
Two situations are considered and compared: (1) wave propagating through the array of cylinders, 
imitating a traditional experimental setup, and (2) wave transmitted from a source located inside the 
ensemble. We show that waves can be blocked from propagation by disorders in the first scenario, 
but the inhibition does not necessarily imply wave localization. Furthermore, the results reveal the 
phenomenon of wave localization in a range of frequencies. 



PACS numbers: 43.20., 71.55J, 03.40K 

When propagating through media containing many 
scatterers, waves will be repeatedly scattered by each 
scatterer, forming a multiple scattering process [Q. Mul- 
tiple scattering of waves is responsible for many fasci- 
nating phenomena such as random laser 0], electronic 
transport in impure solids ||, and photonic or acoustic 
bandgaps Under proper conditions, multiple scat- 

tering leads to the unusual phenomenon of wave local- 
ization. That is, waves in a randomly scattering medium 
are trapped in space and will remain confined around the 
initial transmitting site until dissipated. 

Over the past twenty years, tremendous efforts have 
been devoted to the investigation of the localization 
phenomenon of classical waves in random media (e. g. 
Rcf. J7|-pd|). Observation of classical wave localization 
is a difficult task, partially because suitable systems are 
hard to find and partially because observation is often 
complicated by such effects as absorption and attenua- 
tion. In most previous experimental studies, the appa- 
ratus is set up in such a way that waves are transmitted 
at one end of a scattering ensemble, then the scattered 
waves are recorded either on the other end to measure the 
transmission or at the transmitting site to measure the 
reflection from the sample. The results are subsequently 
compared with the previous theory to infer possible lo- 
calization effects. In this way, observations of wave local- 
ization effects have been reported, for example, for mi- 
crowaves jj), acoustic waves fuJl , and arguably for light 
[fl2| |i~4f , showing propagation inhibition and an exponen- 
tial decay in wave transmission along the propagation 
path, the indications for wave localization, and tending 
to support the prevailing view that waves are localized in 
two dimensional (2D) systems with any amount of disor- 
ders for all frequencies ]l5| ]. 

The purpose of this Letter is twofold. First, we would 
like to point out that traditional experimental methods 
have uncertainties in discerning localization effects, as 
the observation can be obscured by effects like reflection 
and deflection. These effects attenuate waves, resulting 
in a similar decay in transmission and thus making the 
data interpretation ambiguous. We show that while wave 



localization does lead to an inhibition in wave propaga- 
tion, but the propagation inhibition does not necessarily 
imply wave localization. In other words, it is necessary 
to differentiate the situation that waves are blocked from 
transmission from the situation that waves can be actu- 
ally localized in the medium. Second, we show that in 
the system we study, although waves are not localized 
for all frequencies, wave localization is evident in a range 
of frequencies and when there is a sufficient density of 
random scatterers. 

The model in this Letter is acoustic propagation in wa- 
ter containing many parallel air-filled cylinders. Different 
from the common approach that derives approximately a 
diffusion equation for the ensemble-averaged energy, our 
method is to solve the wave propagation from the funda- 
mental wave equation, without resort to approximations. 
The model has been studied previously for the coherent 
behavior of acoustic propagation Jl6| and the acoustic 
complete bandgaps JlTj ] . 

Consider N straight cylinders located at n with i = 
1, 2, • ■ • , N to form either a random or regular array. An 
acoustic line source transmitting monochromatic waves 
is placed at r s . The scattered wave from each cylinder 
is a response to the total incident wave composed of the 
direct wave from the source and the multiply scattered 
waves from other cylinders. The final wave reaching a 
receiver located at f r is the sum of the direct wave from 
the source and the scattered waves from all the cylinders. 
Such a scattering problem can be solved exactly, following 
Twersky |l8|]. While the details are in [fl9|| , the essential 
procedures are summarized below. 

The scattered wave from the j-th cylinder can be writ- 
ten as 
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where k is the wavenumber of the medium, is the 

n-th order Hankel function of the first kind, and oV-r., 
is the azimuthal angle of the vector r — fi relative to the 
positive x-axis. The total incident wave around the i-th 
cylinder (i = 1, 2, • • • , JV; i ^ j) is the summation of the 
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direct incident wave from the source and the scattered 
waves from all other scatterers, and can be expressed as 
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To solve for A % n and B l nl we express the scattered wave 
p s (r, r"j), for each j ^ i, in terms of the modes with re- 
spect to the i-th scatterer by the addition theorem for 
Bessel functions p(| . The resulting formula for the scat- 
tered wave p s (r, fj) is 



p s (r,f-)= J2 
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with 
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The direct incident wave around the location of the i-th 
cylinder can be expressed in a Bessel function expansion 
with respect to coordinates centered at fj 
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with the known coefficients 57 



Using equations (||) , (||) and , we have 
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The boundary conditions state that the pressure and the 
normal velocity be continuous across the interface be- 
tween a scatterer and the surrounding medium, leading 
to 



(7) 



where T l n are the transfer matrices relating the acoustic 
properties of the scatterers and the surrounding medium, 
and have been given by Eq. (21) in Ref. fl9|| . 

The coefficients AL and B 3 n can be inverted from 
Eqs. (H|), (||), and (Q). Once A l n are determined, the 
transmitted wave at any spatial point is given by 
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The acoustic intensity is represented by the squared mod- 
ule of the transmitted wave. 



In the following computation, we assume N uniform 
air-cylinders of radius a. The fraction of area occupied 
by the cylinders per unit area is (3. The average distance 
between nearest neighbors is therefore d = (7r//3) 1 / 2 a, 
which is also the lattice constant for the corresponding 
square lattice array. Two situations are considered: (1) 
wave propagating through the array of cylinders, labelled 
hereafter as the 'Outside' situation that imitates the tra- 
ditional experimental setup, and (2) wave transmitted 
from a source located inside the ensemble, labelled here- 
after as the 'Inside' situation. Both cases are illustrated 
by Fig. [I]. For the 'Outside' case, all cylinders are ran- 
domly placed within a rectangular area with length L 
and width W . The transmitter and receiver are located 
at some distance from the two opposite sides of the scat- 
tering area. For the 'Inside' situation, all cylinders are 
placed within a circle of radius L with the transmitting 
source located at the center and the receiver located out- 
side the scattering cloud. In the computation, the acous- 
tic intensity is normalized in such a way that its value 
equals unity when there are no scatterers present; thus 
the uninteresting geometrical spreading effect is naturally 
eliminated. We scale all lengths by the parameter d, and 
the frequency in terms of non-dimensional ka; in this way, 
the computation becomes non-dimensional. 
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1. Conceptual layout: (a) acoustic propagation 
through a cloud of cylinders; (b) acoustic transmission from 
a line source located inside a cylinder cloud. 

A set of numerical computations has been performed 
for various area fractions (3, numbers N, and dimension- 
less frequency ka. The major controlling parameter is [3. 
The transmitter and receiver are placed at a distance of 
2d from the sample; in fact, we found that as long as we 
keep the symmetry the results remain qualitatively un- 
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changed as the positions of transmitter and receiver vary. 
Though it is the line source that is used in computation, 
the features hold even when a beamed cylindrical plane 
wave is used. 
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FIG. 2. The middle and right panels, referring to the 
'Inside' and 'Outside' cases respectively, show the normal- 
ized acoustic transmission versus frequency in terms of 
non-dimensional ka. Here the comparison is made between 
the results from the corresponding square arrays (solid lines) 
and from the complete random array of cylinders (dotted 
lines). Left panel: The band structures computed by the 
plane wave expansion method. 

Fig. H presents typical results for the transmitted in- 
tensity as a function of ka for the two situations in Fig. [l], 
with (3 = 1CP 3 and N = 200. For reference, we also plot 
the band structure, following [p|Ji^], and the transmis- 
sion for each corresponding square lattice with the same 
(3. Here is shown that for both situations, there is a sig- 
nificant transmission reduction regime from ka = 0.007 
to about 0.022, which is roughly coincident with the 
complete bandgap shown by Fig. ^](a), and within this 
regime, the transmission is less inhibited compared to 
that from the corresponding square lattice arrays. For 
the 'Inside' case, the reduction regime for the random 
scattering is identified as the localization range. This re- 
duction regime will be widened as (3 increases, but will 
disappear when (3 drops below about 10 -5 . Comparing 
the results from the random arrays and that from the 
corresponding square lattice arrays, a significant differ- 
ence is apparent: outside the severe reduction regime, 
the transmission is reduced by the randomization in the 
'Outside' case, for example at ka — 0.05, but stays nearly 
unchanged in the 'Inside' scenario. More explicitly, the 
randomness tends to block the wave propagation outside 
the gap regime. Traditionally, such propagation inhi- 
bition caused by randomness has been regarded as the 
indication of wave localization. In what follows, we show 
that waves are actually not localized at these frequencies. 

We consider two frequencies as an example: ka = 0.01 
and 0.05. Fig. || presents the results for the random 
ensemble averaged transmission and its fluctuation as a 
function of the sample size at the two frequencies for the 
fixed (3 = 10 -3 . For the 'Outside' case, the width of the 
sample is fixed at W — 20. The sample size is varied by 



adjusting the number of the cylinders. A few important 
features are discovered. 

For ka — 0.01, the transmission decays exponentially 
with the sample size for both 'Inside' and 'Outside' situ- 
ations. In the 'Outside' case, there are two decay slops, 
i. e. -0.0427 and -0.0019. There is a transition regime 
separating the two slops. The transmission starts to de- 
cay with slop —0.0427, followed by a milder decay of slop 
-0.0019. The slop of —0.0427 is nearly the same as the de- 
cay slope of —0.0460 in the 'Inside' case. Obviously, this 
is because at ka = 0.01, waves are localized. This expo- 
nential decay actually indicates that waves are trapped 
or localized near the transmitting source; this is clearly 
shown by the top portion of Fig. |3|(b) . The second slop in 
the 'Outside' case for ka = 0.01 is due to the finite width 
of the sample, as will be discussed later. Inside the lo- 
calization regime, the transmission fluctuation is small, 
as expected from an earlier work jl6). Here we see that 
within the localization regime, wave localization can be 
observed in both 'Outside' and 'Inside' scenarios. 
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FIG. 3. Normalized acoustic transmission and its fluctua- 
tion versus the sample size for two frequencies: the left and 
right panels refer to the 'Inside' and 'Outside' cases respec- 
tively. The estimated slops for the transmission are indicated 
in the figure. Here T = p/po- 

For ka = 0.05, the 'Inside' and 'Outside' scenarios dif- 
fer significantly. While for the 'Outside' case the trans- 
mission decreases exponentially with a slop of 0.0027 
along the path, the transmission in the 'Inside' situa- 
tion does not decrease. For the 'Outside' case, the trans- 
mission fluctuation increases then drops along the path, 
emulating the localization effect. For the 'Inside' case, 
however, the transmission fluctuation representing the 
diffusive intensity increases as more and more scat- 
tering occurs along the path, fully complying with the 
well-known non-localized Milne diffusion. The large fluc- 
tuation implies that the transmission is sensitive to the 
distribution of the cylinders, another indication of the 
non-localization property In fact, the apparent de- 
cay in the 'Outside' case is due to the scattering atten- 
uation that the waves are reflected and scattered to the 
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sides. These results suggest that waves are actually not 
localized for ka = 0.05, and it would be a mistake to 
interpret the exponential decay shown in the 'Outside' 
situation as the indication of wave localization. 

Now we consider the width effect in the 'Outside' sit- 
uation. In Fig. [|, the ensemble averaged transmission 
is plotted as a function of the sample size (L) for three 
widths (W) at two frequencies. 
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ka = 0.05, outside the localization regime, the transmis- 
sion decays nearly exponentially along the path. As the 
width increases, the slop will become smaller, and will 
be saturated to a value of -0.0008. These results indicate 
that in the 'Outside scenario', the path-dependent trans- 
mission behaves similarly for frequencies either inside or 
outside the localization regime as long as the width is 
sufficiently large; therefore the phenomenon of wave lo- 
calization cannot be isolated in this scenario. 

Finally, we stress that whether waves are localized or 
extended is an intrinsic property of the system that is 
supposed to be infinite. This property does not depend 
on the source, and should not depend on boundaries ei- 
ther. We believe that while the source is placed inside 
the medium with increasing sizes, the infinite system 
might be mimicked and the localization property could 
be probed without ambiguity. 
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FIG. 4. Normalized acoustic transmission across a scatter- 
ing array as a function of the sample length for different sam- 
ple widths. The solid line in (b) is the numerical extrapolation 
as the width approaches infinity. 

For frequencies inside the localization regime, the sec- 
ond slop is due to the finite width, and can be interpreted 
as follows. At ka = 0.01, the localization effect is domi- 
nant, thus leading to a rapid decay in transmission along 
the path, giving rise to the first slop. As expected, this 
slop is almost independent of the width. As the sam- 
ple size increases, the directly transmitted waves are al- 
most completely blocked. But, due to the finite width, 
a small amount of the deflected waves around the sam- 
ple sides can still reach the receiver. When increasing 
the sample size (L) for a fixed width (W), this effect 
gradually diminishes along the path, yielding the second 
slop. Increasing width (W) will reduce this finite- width 
effect, and thus the amount of deflected waves reaching 
the receiver will also decrease, as shown. The width ef- 
fect disappears when the width approaches infinity. For 
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